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ATTITUDE CONTROL OF A SPACECRAFT WITH A STRAPDOWN 
INERTIAL REFERENCE SYSTEM AND ONBOARD COMPUTER 


by 

John Hrastar 

Goddard Space Flight Center 


INTRODUCTION 

In recent years, the microminiaturization of electronic components has made possible the manu- 
facture and use of very small, general purpose (GP) digital computers, which can be included aboard 
spacecraft. One of the primary uses for such a computer is data processing. Another important use is 
in conjunction with the attitude control system of the spacecraft. The advent of the small computer 
has made it possible to consider relatively complex control laws for the spacecraft attitude control 
system. All the effects of dynamic coupling and large-angle nonlinearities may be considered at the 
start of the design procedure. Small-angle limits need not restrict the control system design. Thus an 
integrated, three-axis control system may be designed in place of three small-angle single-axis systems. 

Although a great deal of work remains in this area, Mortensen (Reference 1) and Meyer (Refer- 
ence 2) have proposed three-axis control laws. 

Meyer defines the direction cosine matrix relating the spacecraft axes to a fixed inertial reference 
as the output of the attitude control system. This matrix is multiplied by a reference matrix, and the 
control law is formulated from the product or error matrix. The control law he assumes is shown to 
be asymptotically stable in the large. The attitude matrix is assumed to be known at all times. It may 
be determined by some type of inertial sensor, e.g., star trackers or gyros. The matrix formulation and 
multiplication and the control law formulation require a computer. 

In an earlier work, using the Cayley-Rodrigues parameters, Mortensen proposed a control law un- 
der the assumption of a particular formulation for the spacecraft kinematics. These are also known as 
Euler- Rodrigues parameters (by Roberson) and the Gibbs vector components (Reference 3). Mortensen 
assumes the parameters defining the body orientation to be known and shows the control using these 
parameters to be asymptotically stable in the large. The parameters are not easily measured but may 
be determined by solving the first-order differential equations in the parameters and body rates. This 
method is similar to Meyer’s method except that the attitude matrix is never formed explicitly. The 
control law is formed directly from the kinematic parameters. The computation of the kinematic 
parameters and the control law formulation require a computer. 
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Control laws of these types have three characteristics that should prove very valuable in future 
spacecraft: 

(1) Stability for any attitude. Thus, the desired attitude of a spacecraft is unrestricted and may 
be determined from considerations other than stability. 

(2) Three-axis reorientation capability. Slewing all three axes simultaneously results in a rapid 
reorientation. 

(3) Simplicity. Since a three-axis reorientation is simply an extension of a hold or point opera- 
tion, mode switching may be simplified. 

This report is a result of a conceptual study for an advanced Orbiting Astronomical Observatory 
(OAO) spacecraft. The spacecraft considered was similar to the present OAO but included an ex- 
tended telescope. The basic assumptions underlying the study with respect to the attitude control 
system were the presence of a GP, onboard computer and a very precise, strapdown inertial reference 
unit (IRU). The primary purpose of the study was to determine the feasibility of a three-axis control 
system’s pointing and reorienting the spacecraft. Use of the computer and IRU was assumed, and 
implementation of the type of control law proposed by Mortensen was considered. Although an ad- 
vanced OAO is the specific spacecraft considered, the principles may be applied to other spacecraft 
that require a pointing and slew capability. 

The control law is formulated along with the algorithm used to generate this control law within 
the onboard computer. The input to the onboard computer is the pulse output of the strapdown IRU. 
The spacecraft dynamics and the onboard computer were simulated on a digital computer to check the 
effect of sampling interval, slew angle, quantization level, and other variables. These simulations gave 
good results for the parameters of this study. A computer and IRU presently under development are 
acceptable for the type of spacecraft studied. 

The simulation was also used to check the dynamic response of the spacecraft to a reorientation 
maneuver. The spacecraft responded well. If time is an important factor, this type of control law is 
superior to the more conventional systems. 

The in-flight reprograming capability of the onboard computer gives this system a great deal of 
flexibility. This is evidenced by the ability of the spacecraft to slew about an arbitrary axis through 
the center of mass as discussed in a later section. 

It is concluded that this type of system is feasible for an observatory type of spacecraft and offers 
a number of advantages over conventional systems. 

CONTROL LAW FORMULATION 

It is necessary to know the spacecraft attitude at all times to effect a three-axis reorientation. 
Knowledge of the orthogonal transformation matrix relating the spacecraft body axes to a known, 
fixed inertial frame is therefore required. Explicit determination of the entire matrix is not necessary 
as long as the three independent parameters of the kinematic representation are known. 
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One method of parameterizing this matrix is by the Cayley-Rodrigues parameters a, 0, and 7 . See 
Equations A 17, Appendix A. In terms of these three independent parameters, the matrix* (Reference 
3) is 


1 , + a 2 + 0 2 + 7 2 


1 + a 2 - 0 2 - j 2 
2 (a 0 -y) 
2(07 + 0) 


2(a0 + 7 ) 

1 - a 2 + 0 2 - 7 2 
2(07 - a) 


2(a7 - 0 ) 

2(07 + a) 

1 - a 2 - 0 2 + 7 2 


( 1 ) 


Continuous knowledge of a, 0, and 7 determines body attitude at all times with respect to the 
inertial frame. The parameters are not easily directly measured but may be determined by continu- 
ously solving the kinematic differential equation relating a, 0, and 7 to the body rates oj , gj , and gj 

x y z 


a 

1 

0 

ii 

to |» 

7 



1 +a 2 
a0 + 7 
O7~0 


a0~7 
1 + 0 2 
07 + a 


07+0 


CO 

X 

07 ~a 


% 

1 +7 2 




( 2 ) 


Using these parameters, Mortensen has proposed the following control law (Reference 1): 


u x = k 2i w x + + o 2 + 0 2 + 7 2 )a , 

u=k 22 co+k 1 (l+a 2 + 0 2 + 7 2 )0, (3) 

and 

u z = k 23^z + k x (l + a 2 + 0 2 + 7 2 )7 , 


where u x , u y , and u z are control torques and k v k 2V k 22 , and k 23 are constant gains. He shows the 
system with this control law to be asymptotically stable in the large with respect to the origin (a = 0 = 
7 = 0). The dynamic and control equations are developed in Appendix A. This control law is analo- 
gous to a conventional law in that the torque is a function (although nonlinear) of rate and position. 
For small angles (i.e., the spacecraft body axes near the desired attitude), the following approxima- 
tions (Reference 3) hold: 


and 



(4) 


where d v 6 2 , and 0 3 are Euler angles. 

*Capital letters with double overbars are square matrices; those with single overbars are column matrices. 


3 



Therefore, for small angles, the control reduces to a simple linear rate plus position control. 


and 


with /c 3 = kj2. 






22 


« A: 


23 


+ fc 3 0 2 , 
+ k 3 e 3 , 


(5) 


A control law of this type, globally stable for any attitude, includes the ability to point and slew 
without the necessity for mode switching. The inertial reference is defined as a = & = y = 0. When a 
new reference is desired, e.g., pointing at a different star, the transformation relating the present posi- 
tion to the new position is determined. The corresponding parameters a, j3, and 7 are determined and 
sent up to the spacecraft. The spacecraft, using this control law, acts to null a, 0, and 7 and thereby 
moves to the new attitude. This is a reorientation or three-axis slew. 


By means of the strapdown IRU and GP computer, a, j3, and 7 may be continuously updated to 
allow use of this control law. (See Figure 1.) 


ATTITUDE REFERENCE ALGORITHM 

The sensor to be used for operational control of the spacecraft during a reorientation is an IRU 
composed of three inertial gyros. This is a strapdown system with the gyros operating in a pulse rebal- 
ance loop. They are continuously nulled by a series of pulses the size and rate of which provide 


DISTURBANCES 



Figure 1— Control system block diagram. 
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position and rate information. References 4 and 5 describe this type of gyro. The information from 
the IRU must be used in conjunction with the computer to update the parameters a, j3, and 7 through- 
out a reorientation. This updating may be done by solving Equation 2. 

The body rates are not known accurately enough for Equation 2 to be solved in its present form. 
Instead an integration algorithm based on an incremental angle may be used (Reference 5). 

A Taylor series expansion may be used to update a, 0, and 7. The value at time t + h may be de- 
termined from the value at time t and the information received from the gyro during the sampling 
interval h. Only the parameter a(t ) will be considered here, although similar expressions exist for j3(r) 
and 7 (t). 


a(t + h) = a(t) + d{t)h + \ fd(t)h 2 + ^rd{t)h z + . . . . (6) 

2 o 

If only a first-order expansion is used and if a substitution for a(t) is made from Equation 2, 

[ 1 + a\t) ]A x [a( 00(0 - 7(01 A y [a(07(0 + 0(01 A z 

a(t + h) = a(t) + - + — + , (7) 

where A ; = c o t {t)h, i = x, y, and z. The A. terms represent the angular output of the ith gyro during 
the interval h, i.e., a number of pulses of weight q (q is the gyro quantization level, 2.4" of arc in the 
IRU slew mode). The updating therefore may be done on an incremental angle basis. If the second- 
order terms of the Taylor series are kept, terms on the order of h 2 appear. Appendix B lists the 
second-order expansions for all three parameters. 

Updating the parameters in this manner allows the control law to be updated at each sampling 
point also. The body rates used in the control law are simply A ^h. 

A study was conducted to determine the feasibility of this type of algorithm. A digital computer 
simulation of the system was used to help evaluate the errors. The attitude error was defined as the 
difference between the true spacecraft attitude and the attitude computed by means of the Taylor 
expansion. The spacecraft dynamics and true attitude were computed by means of a four-point 
Runge-Kutta routine. The integration step was kept much smaller than the sampling interval h. 

The results of the study are as follows: 

(1) The error in the update algorithm increases with an increase in the sampling interval h. This 
is not a serious constraint because even relatively long sampling intervals gave good results. With a 
second-order Taylor expansion and a 5-s sampling interval, the errors at the end of a 60° slew were on 
the order of 14" of arc. 

(2) The error in the update algorithm increases with the magnitude of the slew angle for the 
longer sampling intervals (2 s and greater). For a sampling interval of 100 ms, the error is apparently 
independent of the magnitude of the slew angle even for very long slews (165°). It remains on the 
order of 2.4" of arc, which is the quantization level of the gyros. For a sampling interval of 1 s, the 
error for a 165° slew is only about 5.3" of arc with a second-order Taylor expansion. 
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(3) For moderate sampling intervals and slew angles (up to 2 s and 90°), the attitude error is 
primarily a function of the gyro quantization level and is on this order of magnitude. This also holds 
true for shorter slews (30°) with longer sampling intervals (5 s). 

(4) The computer round-off error may be kept negligible by the use of a 36 bit, double-precision 
word. The computer considered in the study (see the following section) has the ability to operate on 
either a single- (18 bits) or double-precision basis. 

The estimated times for the update computations are shown in Table 1. These times are based on 
the computer described in the following section and are compatible with the sampling intervals inves- 
tigated. 

These results indicate that use of this type of algorithm for updating the spacecraft attitude 
parameters is feasible. 


Table 1— Estimated time for the update computation. 


Taylor series expansion 

Single-precision 

(ms) 

Double-precision 

(ms) 

First-order 

2 

7 

Second-order 

4 

14 


COMPUTER 

A candidate for the computer required for attitude updating and control law formulation is the 
onboard processor being developed at GSFC.* This computer is being developed primarily as a data 
processor but appears to have the capacity and capability for handling the control equations. Some of 
the features of this computer are 

(1) Postlaunch reprogramming capability. The control laws may be changed in flight if desired. 

(2) Low power. The central processing unit requires 5 W. 

(3) Add time of 7.5 jus. 

(4) Maximum use (80%) of integrated circuits. 

INERTIAL REFERENCE UNIT 

The IRU being developed by MIT for the OAO will apparently meet the attitude reference re- 
quirements. Three operational modes are available: hold, slew, and I-stab. The primary concern of 
this report, three-axis reorientation, requires use of the slew mode. For this study the processing being 


’’Taylor, T., et. al., “A General Purpose On Board Processor for Scientific Spacecraft,” NASA/Goddard Space Flight Center 
Document X-562-67-202, July 1967. 


6 



developed for the early models of the IRU was not considered. The only concern was the basic per- 
formance of the gyros. 

Some of the important characteristics of the IRU in its slew mode are a resolution of 2.4" of arc / 
pulse, a maximum rate of 480" of arc/s, a pulse rate of 200 pps, a compensated drift of 20" of arc / 
orbit ( 1 2" of arc/hr), and a torquing error less than 0.01%. 

NORMAL REORIENTATION 

The system dynamic, kinematic, and control equations were simulated on a digital computer. 

The primary objective of the simulation was to determine performance with respect to a slew or re- 
orientation. The equations reduce to approximately linear equations for small angles (hold or pointing). 
Thus it was not necessary to consider the hold or pointing modes of operation. 

The angle 0 (Appendix A) is a good scalar representation of error, i.e., the difference between the 
actual body orientation and the desired orientation (Reference 2). This is true because a, j3 , and y are 
all zero if and only if 0 is zero. Therefore, 0 was chosen as the primary performance index. For a 
single-axis slew, 0 is equivalent to the Euler angle around the slew axis. 

The only requirements for stability are that the position gain k and rate gains k , k, and be 

P x y z 

positive and nonzero. There are no other requirements with respect to magnitude or linearity. There- 
fore, once the stability criteria are met, the designer may set the gains to meet other criteria such as 
speed of response, damping, etc. 

Response 

Probably the most important criterion for a normal reorientation is the time. The shorter the 
reorientation time, the longer the experimenting time. Another criterion, maintenance of a stable slew 
axis, is discussed in a following section. The time criterion was used in making several computer runs. 

The time for a slew was defined as the interval which extended from the time when initial condi- 
tions were read into the onboard computer to the time when the norm reached its minimum value. 

The norm, defined as 

NORM = + + cuj) + + 0 2 ) 1/2 , (8) 

was usually assigned a minimum of 10~ 4 . Since the norm is always positive, the rates were required to 
be below 10“ 4 rad/s and 0 was required to be below ICf 4 rad (0.34' of arc) before a slew was consid- 
ered complete. Therefore, slew or reorientation time as defined here includes settling time. 

Since a slew is basically a momentum exchange, the gains were set to result in a maximum wheel 
speed during the slew. The rate gains were set to provide good damping at the end of the slew. Both 
type A and type B wheels were used. (See Figure 2.) The characteristics were a stall torque of 0.27 
N-m (0.20 ft-lb), a maximum momentum of 13.6 N-m-s (10.0 ft-lb-s), and a time constant T of 50 s. 
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Figure 2— Reaction wheel torque speed curve. 


(a) type A wheel and (b) type B wheel. 


To compare a three-axis reorientation with the more conventional method using three single-axis 
slews, a number of runs were made with the slew axis eigenvector aligned with a control axis. (See 
Equations A 17, Appendix A.) This is equivalent to a conventional single-axis slew. 

Figure 3 shows the response for four single-axis and two multiaxis slews [eigenvector (0.5, 0.5, 
0.707)]. The maximum slew rate for the single-axis slew is about 0.25 X 10 -2 rad/s. This is deter- 
mined by the 13.6 N-m-s maximum wheel momentum and the 5420 kg-m 2 body inertia (the same for 
all axes). The maximum slew rate for the multiaxis slew occurs when all three wheels are saturated. It 
is about 0.43 X 1(T 2 rad/s. Thus all three wheels are being used to reduce 0. The single-axis and 
multiaxis slews in Figure 3 are not directly comparable because the final orientations are different. A 
comparison will be made later in the report. 

Figure 4 is a plot of total slew time for slews of various lengths. The curves all converge to about 
160 s for a 2° slew. This is apparently due to the end effects, i.e., acceleration at the start and settling 
at the completion of the slew. The slope of the linear portion of the single-axis curves is (0.25 X 10 -2 ) 1 
s/rad. This is determined by the maximum wheel momentum. All the points checked for multiaxis 
slews were below the time required for equivalent length single-axis slews. 

Three points were checked both ways to compare the time required for reorientation by means 
of a multiaxis slew with the time required for the same reorientation by means of three single-axis 
slews. See points D, E, and F on Figure 4. Table 2 summarizes the results. 

In the table, T x is the total time for three single-axis slews. This time is determined by using the 
Euler angles and single-axis slew curves of Figure 4. The time to complete the same reorientation using 
a multiaxis slew is Ty The ratio shows the multiaxis reorientation to be more than twice as fast 

as the conventional slew. This should be expected for at least two reasons: The end effects, which are 
independent of the length of the slew, occur three times for three slews and only once for a single 
slew, and all the wheels contribute simultaneously during a multiaxis slew. This contribution causes 
the total rate by which the error is being reduced to be higher than for a single-axis slew. 

The slew time using a type B motor was slightly longer in each case due to the reduced torque 
available for acceleration and deceleration. The slopes of the linear portions are the same since the 
maximum wheel speeds are the same. 
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Accuracy 

In a hold or point mode the accuracy is primarily a function of the sensors. When the spacecraft 
is holding on the experiment telescope or roll axis star tracker, there is no problem. When it is holding 
on the IRU, the primary error will be due to the gyro drift. This is expected to be about 1 2" of arc/hr. 
The desired accuracy determines the update frequency when the spacecraft is holding on the IRU. 
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Figure 4-Slew time. 


Table 2— Reorientation times. 



Euler angles 
(rad, deg) 

Eigenvector comp. 

0 

*1 

T 
1 3 




*2 

^3 

c. 



(rad, deg) 

(s) 

(S) 


D 

0.1745, 

10 

0.2745, 

15.7 

0.1745, 

10 

0.5187 

0.6797 

0.5187 

0.3797, 

21.8 

700 

260 

2.69 

E 

0.523, 

30 

0.523, 

30 

0.523, 

30 

0.6546 

0.3782 

0.6546 

0.9722, 

55.7 

1080 

420 

2.57 

F 

1.045, 

60 

1.045, 

60 

1.045, 

60 

0.6946 

0.1869 

0.6946 

2.0327, 

116.5 

1680 

790 

. _ J 

2.13 


There are three sources of error present during a reorientation. They are the attitude updating, 
gyro drift, and the torquing inaccuracy of the gyros. For a sampling interval of 100 ms, the attitude 
update error may be kept on the order of the gyro quantization level, which is 2.4" of arc. The reori- 
entation error due to gyro drift will be small because most slews will take less than 20 min to complete 
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The compensated drift during this time is only 3 to 4" of arc. There is an error accumulation due to 
the torque rebalance loop when the gyro is being torqued. The torque rebalance pulses will contain 
some errors in size and shape. This will cause an attitude error buildup over a long slew. The torquing 
inaccuracy is estimated to be 0.01%. This represents an error of 0.0165° (1' of arc) over a 165° slew. 
Therefore, this is the major attitude error in a large angle reorientation. 

The accuracy of the system need not be as high as that required for fine pointing. The objective 
of the reorientation is the acquisition of a new experiment star. Therefore, the primary factor in the 
determination of the required slew accuracy is the field of view of the experiment telescope or the roll 
axis star tracker. If the total reorientation error is well within the field of view of the sensor, it is fea- 
sible to reorient the spacecraft by the method described. With an experiment that has a field of view 
±4' of arc, it does appear feasible to move from fine pointing to fine pointing by means of a single 
large-angle three-axis reorientation. 

Disturbances 

The effect of environmental disturbance torques during a slew is to increase the total momentum 
of the system. The reorientation trajectory with disturbances will in general be different from the 
undisturbed trajectory. However, the final reference as stored in the computer and IRU is not affected; 
the system will null at the same reference as the undisturbed system. The time for the reorientation 
will depend on the magnitude and direction of the disturbances and may be shorter or longer than that 
required for an undisturbed reorientation. Although an extensive torque disturbance analysis was not 
conducted, some computer runs were made with disturbance torques. No gross changes in system per- 
formance were noted. The additional time required for one 165° slew was on the order of 1 min with 
disturbances on the order of 0.0001 to 0.0005 N-m (1000 to 5000 dyn-cm). 

FIXED-AXIS REORIENTATION 

When the spacecraft is reoriented from one inertial attitude to another, as discussed in the pre- 
vious section, the primary criteria are time and accuracy. Once the attitude parameters are loaded in 
the computer, the spacecraft follows a trajectory determined in part by the controller gains to the new 
attitude. In this case, the instantaneous axis of rotation (the +1 eigenvector of the attitude reference 
matrix) is not inertially fixed but moves throughout the reorientation time. 

There are times when it is important to maintain rotation about a fixed, although arbitrary, axis. 
This would be the case in an acquisition mode in which it is necessary to slew about an axis so that a 
particular sensor will pick up a particular star. The axis of rotation must remain inertially fixed; i.e., 
the rotation eigenvector must remain fixed, to assure the sensor will eventually have the star in its field 
of view. This may be accomplished by commanding a reorientation around this eigenvector and limit- 
ing the motor torques and maximum momenta by the ratio of the eigenvector components. The body 
rates are then proportional to the eigenvector components, and the total angular velocity vector re- 
mains collinear with the slew axis (eigenvector) throughout the entire slew. 

A computer simulation confirmed the stability of the eigenvector position for this type of 
reorientation. 


11 



DISCUSSION 

This study was not a general study of various control laws and/or kinematic representations. The 
choice of Mortensen’s control law with the Euler-Rodrigues parameters appeared to be natural for this 
problem. Further study may show it is advantageous to use Meyer’s method when the entire transfor- 
mation matrix is updated. Basically there is little difference between the two approaches. It may be 
shown that the control laws are very similar. 

These types of control laws may be used with other types of sensors. In particular, the attitude 
matrix defined by Meyer as the position output of the system is independent of the sensors used to 
describe its elements. Star trackers, as well as gyros, may be used. However, it appears that a strap- 
down IRU as considered here offers the most advantages for large-angle reorientations. The use of 
gimbaled star trackers has two disadvantages. Due to the limited gimbal rotation, large-angle reorienta- 
tions cannot be completed without switching trackers during the slew. This introduces additional con- 
straints on the reorientation maneuver. Secondly, the attitude matrix is more complex because of the 
presence of the trigonometric functions resulting from the tracker gimbals. 

Because of the gyro drift it seems necessary for the spacecraft to hold its position during the ex- 
periment time by means of optical trackers, either the experiment telescope or the roll axis star tracker. 
Holding the spacecraft position with an optical tracker requires a relatively simple control system and 
does not require a computer. Therefore, during this time the computer may be free to process the 
experiment data. During a reorientation when no experiment is being carried out, the computer may 
be used completely for the control system. Although the computer-controlled system may easily han- 
dle both the hold and reorientation mode, it may be more efficient to use it in this manner. 

The in-flight programing capability of the computer may be used to advantage. The control gains 
may be set initially to carry out a constant eigenvector slew for initial acquisitions. Once this is com- 
plete, the gains may be changed to allow a more rapid slew for normal reorientations. 

CONCLUSIONS 

A computer control system that performs as well during a slew as during pointing is feasible for 
an advanced OAO. With an onboard digital computer and strapdown IRU, the system may complete 
large-angle three-axis reorientations with the same control law used for pointing. The computer speed 
is such that a second-order Taylor series may be used to update the spacecraft attitude during a reori- 
entation. The errors due to the attitude updating may be kept on the order of the gyro quantization 
level. The primary error in a large-angle reorientation is due to the gyro torquing inaccuracy. This 
error, however, is not prohibitive. 

A GP computer now under development, the On Board Processor being developed at GSFC, 
appears to meet the computational requirements. 

An IRU being developed for GSFC by MIT appears capable of meeting the sensor requirements. 
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A three-axis slew is superior to three single-axis slews when time is important. A three-axis slew 
was shown to be more than twice as fast as three single-axis slews for the reorientations that were 
checked. It is also simpler and requires fewer commands at the start of the slew. 

If the proper gain settings are used, a slew may be maintained about any constant, arbitrary axis. 

It is possible to use a conventional linear rate and position control when pointing and to employ 
the computer primarily for reorientation. The computer could then be used for data processing when 
the spacecraft is in a pointing mode. 
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Appendix A 


DYNAMIC AND STABILITY ANALYSIS 


The spacecraft dynamic and control equations are developed in this appendix. Although Mortensen 
did not use reaction wheels for control, the method for proving stability closely parallels his develop- 
ment (Reference 1). 

If the control torque on the spacecraft is assumed to be developed by reaction wheels, the dynamic 
equations (Reference 3) are 

K * V, - + wz - ", H ,„> ‘ ■ 

it, + v, ( A » 

and 

K i-O y H WX ) “ ^2 H wz , 


where H v , H„, and H, are the spacecraft momentum about x, y, and z (principal axes) minus wheel 
momentum ( H f = /.cc., i = x, y, and z); I x , I y , and I z , the spacecraft inertia about x, y, and z; H wx , 
H wy , and H wz , the wheel momentum about x, y, and z; and M x , M y , and M v the external torques 
about x, y, and z. Assume M x = M y = M z = 0. During a slew the external torques are small compared 
to the control torques, and the assumption should be valid. If the equation is put into matrix form. 
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H„ 


H. 


H. 


Tz 


-H 


Ty 


-H. 


Tz 


H, 


Tx 


H. 


Ty 


-H, 


Tx 


0 




H„ 


H. 


H.. 


H, 


wy 


H 


(A2) 


where H Tj = H t + H wj , i = x, y, and z. If Equation A2 is combined with Equation 2, 
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I 



(A3) 


If a control law of the type shown in Equations 3 and a reaction wheel with a torque speed curve as 
shown in Figure 2(a) is assumed, the control equations are 


H wx = k x H x + k p 


1 +q 2 
2 


a + k 





where k v , k , and k are rate gains and k is position gain. 

x y z p 

The total system is ninth order and requires nine state variables for a complete description of the 
trajectory in state space. The state variables are the three body momenta, three parameters specifying 
body position, and the three total momenta. The equilibrium position of the spacecraft can be seen 
to be independent of the system total momenta. Thus there is a three-dimensional subspace every 
point of which is an equilibrium point for the spacecraft. Therefore, for spacecraft stability analysis 
the six-dimensional vector representing the body position and rate may be used in place of the nine- 
dimensional system vector. The total momentum in this formulation may be considered a time-varying 
parameter. 
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This is in the form 


X = F(X)X 


where X - [H x , H y , H z , a, (3, y\ T and F(X) is a nonlinear system matrix. Define the matrices K 
as 


and 


K = diag [~k x , -lc y , -k gt 0, 0, 0] 


R = F- K . 

Choose as a candidate V function the quadratic form 

v= x T dx 


where 


G = diag 


Hi , . . ' 

I »/ ’/ » k p’ k p> k P 
. x y z 


Since G is a constant, it follows that 

v = i T dx + x t gx . 

If Equations A7 are substituted into A6 and then into A9, Equation A9 becomes 
v = x T R T dx + x T f T dx + x T dKx + X T (5EX . 


(A5) 


(A6) 
and R 

(A7) 

(A8) 

(A9) 

(A10) 
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(All) 


Note that GR is skew symmetric and that G = G T . Therefore, 

UK = - (UKf = -K t G t = -K T U . 

This reduces Equation A10 to 

V = X T 1 T UX + X T UK.X . (A 12) 

K T U = K T U T = (UK) t =GK (A 1 3 ) 


is also true. Therefore, 


Thus, if 


V= 2 X t GKX 


% -*'*f*,*rA- 


k x > 0 , k y > 0 , and k 2 > 0 , 


V is negative semidefinite. Since 


V 




H 2 Z 

+ — + k p a? + k p p 2 + k p y 2 


(A 14) 
(A15) 


(A16) 


is positive definite (k > 0), the origin is stable. It may be shown that the origin is not only stable but 
asymptotically stable; i.e., V + 0 as t * °°. 


Since the relations 


and 


0 

a = C x tan - , 

P = Cy tan — , (A17) 

0 

7 ~C Z tan — , 


are valid (Reference 3), where C x , C , and C z are the components of the +1 eigenvector of the orthog- 
onal transformation parameterized by a, (3, and y and 0 is the angle around the +1 eigenvector that the 
body frame is rotated from the inertial frame, Equation A16 may be written as 


H l H l H l J 

v*T + f + T +k p m T 

x y z 


(A18) 


Therefore, for V to be bounded it is necessary that the body angular momentum be bounded and 
that 0 be less than n. The total system momentum must be within the total capacity of the wheels. It 
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will be stored in the wheels at the completion of the reorientation. The latter constraint (<p < 7r) is not 
a practical restriction on the system since it is easily satisfied. With the boundedness of V established, 
it may be concluded that V is a valid Lyapunov function and the stability conclusions are valid. 

Asymptotic stability in the large has not yet been proven in general when a type B reaction wheel 
is used. However, two important special cases may be shown to be asymptotically stable. One is the 
case in which the total system momentum is zero. The other is the linearized system near the origin. 

When a type B reaction wheel is used, the additional back emf terms ~H wi lT m , i = x, y, and z, are 
added to the respective control equations (Equations A4). When the total system momentum is zero, 
H { = ~H wV Thus, the system is reduced to sixth order and is easily proven stable in a manner analo- 
gous to that used in the preceding proof. 

When the system equations are linearized near the origin, the result is three uncoupled linear 
second-order systems. These are easily proven stable by conventional methods. Any stored momen- 
tum represents a shift to a new origin. 
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Appendix B 


PARAMETER UPDATE EQUATIONS 

The equations for updating the parameters a, and 7 by means of a second-order Taylor series 
are 

a n + 1 = a n + |[(l + a Z) A x + (“A ~ %,) A y + (%7 n + P n ) A z ] 

+ \ [ a „ 0 + a n) A x + Pn (“A ~ 7„) A y + 7 n (a„7„ + P„) A ] 

+ (2 afr n +p n ~ a n 7 „)A x A y + (2 a 2 n y n + 7 „ + A x A z 

+ KU + ^ 2 -7„ 2 )A/J’ 

=Pn + y[(“A + 7„)A X + (1 +tf)A, + (0„7„ -a„)Aj 

+ i [ a « ( a A + T„) A 2 + P n (l + ^)Aj + 7 „ (P n 7 n ~ a„) A 2 2 

+ (2a„0 2 + a„ + P n 7 n ) A x A y + (2a„P n y n - a 2 + 7 2 ) A jc A z 

+ ( 2 ^7„+7 w - a A) A y A J ’ 

and 

T„ + i = 7„ + ^[(a„7„ * P n ) A x +(P n 7 n + a n) A y + (l + 7„ 2 )A z ] 

+ i [ a « ( a «%, - Pn) A x + Pn(P n y n + a n) A y + T„ (l + 7„ 2 )A 2 

+ ( 2 a A%* +a l~ fi) A x A y + ( 2a nV 2 n +a n~ A x A y 
+ (2^7 2 +^+ a „7„)A y A z ] . 
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